Decoherence in a dynamical quantum phase transition 
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Motivated by the similarity between adiabatic quantum algorithms and quantum phase transi- 
tions, we study the impact of decoherence on the sweep through a second-order quantum phase 
transition for the prototypical example of the Ising chain in a transverse field and compare it to 
the adiabatic version of Grovers search algorithm, which displays a first order quantum phase tran- 
sition. For site-independent and site-dependent coupling strengths as well as different operator 
couplings, the results show that (in contrast to first-order transitions) the impact of decoherence 
caused by a weak coupling to a rather general environment increases with system size (i.e., number 
of spins/qubits). This might limit the scalability of the corresponding adiabatic quantum algorithm. 

PACS numbers: 03.67.Lx, 03.65. Yz, 75.10.Pq, 64.60.Ht. 



I. INTRODUCTION 



A. Quantum Phase Transitions 



In contrast to thermal phase transitions occurring 
when the strength of the thermal fluctuations reaches 
a certain threshold, during recent years, a different class 
of phase transitions has attracted the attention of physi- 
cists, namely transitions taking place at zero temperature 
An analytic non-thermal control parameter such as 
pressure, magnetic field, or chemical composition is var- 
ied to access the transition point. Despite the analytic 
form of the order parameter, the ground state of a sys- 
tem changes non-analytically. There, order is changed 
solely by quantum fluctuations, hence the name quan- 
tum phase transition (QPT). Let us consider a quantum 
system (at zero temperature) described by the rlamilto- 
nian H depending on some external parameter g. At 
a certain critical value of this parameter g Cl the sys- 
tem is supposed to undergo a phase transition, i.e., the 
ground state \^<(g)} of H(g) for g < g c is strongly dif- 
ferent from the ground state \iff > (g)) of H(g) for g > g c . 
For example, \^<(g)) and \ } i> > (g)) could have different 
global/topological properties (such as magnetization) in 
the thermodynamic limit. Therefore, a quantum phase 
transition can also be defined as a non-analyticity of the 
ground state properties of the system as a function of 
the control parameter. If this singularity arises from a 
simple level crossing in the ground state, see Fig. (JTJ)- 
(a), then we have a first-order quantum phase transi- 
tion. The situation is different for continuous transi- 
tions, where a higher-order discontinuity in the ground 
state energy occurs. Typically, for any finite-size system 
a transition will be rounded into a crossover, this is noth- 
ing but an avoided level-crossing in the ground state, see 
Fig. (JTJ- (fo) . Continuous transitions can usually be char- 
acterized by an order parameter which is a quantity that 
is zero in one phase (the disordered) and non-zero and 
possibly non-unique in the other (the ordered) phase. If 
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Figure 1: (Color Online) Sketch of the lowest eigenvalues of a 
Hamiltonian H(g) as a function of some external parameter g 
for a first order QPT. At the critical point g = g c , the ground 
state changes from \9<(g)) to \^>(g)). (a) A level-crossing. 
(6) An avoided level-crossing. 



the critical point is approached, the spatial correlations 
of the order parameter fluctuations become long-ranged. 
Close to the critical point the correlation length T di- 
verges as T _1 cx A|g — <7 c | a , where is a critical expo- 
nent and A is an inverse length scale of the order of the 
inverse lattice spacing. Let A denote the smallest energy 
excitation gap above the ground state. In most cases, it 
has been found [l| that as g approaches g c , A vanishes as 
A cx T~ z cx A z \g — g c \ 0z , where z is the dynamic critical 
exponent. This poses a scalability problem for adiabatic 
ground state preparation schemes (see below), as these 
require a nonvanishing energy gap A. 



B. Adiabatic Quantum Computation 

Unfortunately, the actual realization of usual sequen- 
tial quantum algorithms (where a sequence of quantum 
gates is applied to some initial quantum state, see, e.g., 
[2|) goes along with the problem that errors accumu- 
late over many operations and the resulting decoherence 
tends to destroy the fragile quantum features needed for 
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the computation. Therefore, adiabatic quantum algo- 
rithms have been suggested |3|, where the solution to a 
problem is encoded in the (unknown) ground state of a 
(known) Hamiltonian. Since there is evidence that, in 
adiabatic quantum computing the ground state is more 
robust against decoherence - the ground state cannot 
decay and phase errors do not play any role, i.e., er- 
rors can only result from excitations 0, H| _ this scheme 
offers fundamental advantages compared to sequential 
quantum algorithms - a sufficiently cold reservoir pro- 
vided. Suppose we have to solve a problem that may 
be reformulated as preparing a quantum system in the 
ground state of a Hamiltonian H( . The adiabatic theo- 
rem then provides a straightforward method to solve 
this problem: Prepare the quantum system in the (known 
and easy-to-prepare) ground state of another Hamilto- 
nian Hq. Apply Hq on the system and slowly modify it to 
Hf. The adiabatic theorem ensures for a non- vanishing 
time-dependent energy gap that if this has been done 
slowly enough, the system will end up in a state close to 
the ground state of Hf . Therefore, a measurement of the 
final state will yield a solution of the problem with high 
probability. 

Furthermore, adiabatic quantum algorithms display 
a remarkable similarity with sweeps through quantum 
phase transitions @, Q ■ For all interesting systems dis- 
cussed later in this article, adiabatic quantum compu- 
tation inherently brings the quantum system near to a 
point which is similar to the critical point in a quantum 
phase transition. As an example for the deformation of 
the Hamiltonian, one can consider the linear interpola- 
tion path between these two Hamiltonians 



H(g) = [l-g(t)]H +g(t)H i , 



(1) 



with <?(0) = and g(T) = 1, where T is the total evolu- 
tion time or the run-time of the algorithm. We prepare 
the ground state of Hq at time t = 0, and then the state 
evolves from t = to T according to the Schrodinger 
equation. At time T, we measure the state. Accord- 
ing to the adiabatic theorem, if there is a nonzero gap 
between the ground state and the first excited state of 
H(g) for all g £ [0,1] then the success probability of the 
algorithm approaches 1 in the limit T — > oo. How large 
T should actually be is roughly given by [32| (for a more 
detailed discussion see, e.g. [34|, | 



max 



se[o,i] 



;i,.< 



0,9) 



lin s6[o,i] [Ei(g) - E (g)Y 



(2) 



where Eg(g) is the lowest eigenvalue of H(g), Ei(g) is 
the second-lowest eigenvalue, and \0,g) and \l,g) are the 
corresponding cigenstates, respectively. Somewhere on 
the way from the simple initial configuration Hq to the 
unknown solution of some problem encoded in Hf , there 
is typically a critical point which bears strong similari- 
ties to a quantum phase transition. At this critical point 
the fundamental gap (which is sufficiently large initially 
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Figure 2: (Color Online) Sketch of the two lowest energy 
eigenvalues of the Grover Hamiltonian (a). In the continuum 
limit, this corresponds to the time evolution of the energy 
landscape for a first-order transition (b). The green dot in 
the energy landscape denotes the ground state. 



and finally) becomes very small, see, e.g., Fig. @. Near 
the position of the minimum gap, the ground state will 
change more drastically than in other time intervals of 
the interpolation. In the continuum limit, one would gen- 
erally expect that the minimum value of the fundamental 
gap in adiabatic computation will vanish identically and 
that the ground state will change non-analytically at the 
critical point. This is similar to what happens in quan- 
tum phase transition when g approaches g c . Based on 
this similarity, it seems @ that adiabatic quantum al- 
gorithms corresponding to second-order quantum phase 
transitions should be advantageous compared to isolated 
avoided level crossings (which are analogous to first-order 
transitions) . A brief review of this idea comes in the fol- 
lowing section. 



II. EXAMPLES 

A. First-Order Transition — Grovers Algorithm 

An adiabatic version of Grovers algorithm [91] is defined 
by the Hamiltonian 



H(g) = (1 - g)H + gH t 



(3) 



where the initial Hamiltonian is given by 
Hq = 1 — | in) (in | with the initial superposition state 
|in) = J2x=o \ x ) I^E and D = 2 N denotes the di- 
mension of the Hilbert space for N qubits. The final 
Hamiltonian reads Hf = 1 — \w) (w\, where \w) denotes 
the marked state. In this case, the commutator is very 
small [Hq,H{] = (| in) (w\ — \w) (in|)/-\AD and one can 
nearly diagonalize both Hamiltonians simultaneously 
and the ^-dependent spectrum will consist of nearly 
straight lines - except near g c = 1/2, where we have an 
avoided level-crossing, see Fig. ©. In the continuum 
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Figure 3: The ground-state energy of the Grover Hamilto- 
nian and its first-order derivative. The discontinuity in the 
first-order derivative of the ground state suggests the first- 
order quantum phase transition in adiabatic quantum search 
algorithm. 
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Figure 4: (Color Online) Sketch of the two lowest energy levels 
of the Ising Hamiltonian (left) given in Eqn. ([5| and the time 
evolution of the energy landscape for a second-order transition 
(right). A symmetry-breaking transition corresponds to the 
deformation of the energy landscape. The green dot in the 
energy landscape denotes the ground state. 



limit of A" — > oo, this corresponds to a first-order 
quantum phase transition from |in) = |— > ■ ■ ■ — >) to 
\w) = |tl " " ■ t44-)i f° r example, at the critical point 
g c = 1/2. Such a first-order transition is characterized 
by an abrupt change of the ground state - |in) for g < g c 
and \w) for g > g c - resulting in a discontinuity of a 
corresponding order parameter, see Fig. ([3]), 

<*<»)lf !*><.» -f • <*> 

In contrast to the conventional order parameters, for lin- 
ear interpolations in Eqn. ([3]) the operator dH/dg treats 
both phases symmetrically. Since [Hi, Hp] =^ (for non- 
trivial systems), dH/dg is off-diagonal in either phase 
and thereby plays an equivalent role such as e.g. magne- 
tization. 

First-order quantum phase transitions are typically as- 
sociated with an energy landscape pictured in Fig. (|2|). 
where the two competing ground states are separated by 
an energy barrier throughout the interpolation. In or- 
der to stay in the ground state, the system has to tun- 
nel through the barrier between the initial ground state 
| in) and the final ground state \w) during the quantum 
phase transition. The natural increase of the strength 
of the barrier with the system size N yields to the tun- 
neling time which scales exponentially with the system 
size. Specifically, the optimal run-time for the adiabatic 
search algorithm behaves as T = O(Vl)) = O (2 N / 2 ) @. 
The observation that this first order QPT is associ- 
ated with an exponentially small energy gap right at 
the avoided crossing can be generalized [l9| to local 
Hamiltonians: The two-dimensional subspace of the 
avoided crossing is spanned by the eigenstates \w<(g)} 
and \w > (g)) that become degenerate at g — g c . Due 
to their macroscopic distinguishability, the overlap be- 
tween these states is exponentially small, which for lo- 



cal Hamiltonians also transfers to the matrix element 
(tu<(0 c )| H(g c ) |iu>(ff c )) = O (exp{-D». Consequently, 
one can conclude from the eigenvalues of H(g) in this 
two-dimensional subspacethat also the minimum energy 
gap will become exponentially small in this case. 

Therefore, the abrupt change of the ground state and 
the energy barrier between the initial and final ground 
states suggest that the first-order transitions are not the 
best choice for the realization of adiabatic quantum al- 
gorithms Q. Thus, it would be relevant to study higher- 
order quantum phase transitions for this purpose. 



B. Second-Order Transition — Ising Model 

The one-dimensional quantum Ising model is one of 
the two paradigmatic examples [l| for second-order quan- 
tum phase transition (the other is Bosc- Hubbard model) . 
Of these two, only the former model is exactly solvable 
0, [l(| - the Ising model in a transverse field is a special 
case of the XY model (which can also be diagonalized 
completely). This model has been employed in the study 
of quantum phase transitions and percolation theory [l| , 
spin glasses p], [TTj] , as well as quantum annealing [l2r - 
Il4| etc. Although its Hamiltonian is quite simple, the 
Ising model is rich enough to display most of the basic 
phenomena near quantum critical points. Furthermore, 
the transverse Ising model can also be used to study the 
order-disorder transitions at zero temperature driven by 
quantum fluctuations [l|, [l2|]. Finally, two-dimensional 
generalizations of the Ising model can be mapped onto 
certain adiabatic quantum algorithms (see, e.g., [la]). 
However, due to the evanescent excitation energies, such 
a phase transition is rather vulnerable to decoherence, 
which must be taken into account (l6| . 

The one-dimensional transverse Ising chain of N spins 
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Figure 5: The ground-state energy for the Ising model and its 
first and second derivatives. The discontinuity in the second- 
order derivative of the ground state suggests a second-order 
quantum phase transition. 



exhibits a time-dependent nearest-neighbor interaction 
g(t) plus transverse field B{t) = 1 — g(t) 

N 

H sys (t) = - X) {[1 " 9(t)] °* + 9(t) ^M+i } . ( 5 ) 
i=i 

where <Tj — (a* , a^j are the spin- 1/2 Pauli matrices 
acting on the jih qubit and periodic boundary condi- 
tions ctn+i = Gi a- re imposed. This Hamiltonian is in- 
variant under a global 180-degree rotation around the 
crj-axes (bit flip) which transforms all qubits according to 
cr| — > -<j|. Choosing g(0) = and g(T) — 1 where T is 
the evolution time, the quantum system evolves from the 
unique paramagnetic state |in) = |— >— >— > . . . ) through a 
second-order quantum phase transition (see figure [5j at 
g CI = 1/2 to the symmetrized combination in the two- 
fold degenerate ferromagnetic subspace, see also Fig. @] 
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(6) 



At the critical point g CI the excitation gap vanishes (in 
the thermodynamic limit N — > oo) and the response time 
diverges. As a result, driving the system through its 
quantum critical point at a finite sweep rate entails inter- 
esting non-equilibrium phenomena such as the creation 
of topological defects, i.e., kinks [171 ] . 

Since the initial ground state |in) reflects the bit- flip in- 
variance of the Hamiltonian ([5]) whereas the final ground 
state subspace |f . . . f) an d |i • • • -V) breaks this symme- 
try, we have a symmetry-breaking quantum phase transi- 
tion. Typically, such a symmetry-breaking change of the 
ground state corresponds to a second-order phase transi- 
tion 0. For such a transition, the ground state changes 
continuously and the energy barrier observed in first- 
order transitions is absent: Initially, there is a unique 
ground state but at the critical point, this ground state 




Figure 6: (Color Online) Sketch of the time evolution of the 
energy landscape for a symmetry-breaking quantum phase 
transition which is of first order. The green dot denotes the 
ground state. 



splits up into two degenerate ground states which are the 
mirror image of each other. Therefore, the ground state 
does not change abruptly in this situation and the sys- 
tem does not need to tunnel through a barrier in order to 
stay in the ground state, see Figs. (JDd)- Consequently, we 
expect that in this closed quantum system should 

find its way from the initial to the final ground state 
easier. This expectation is confirmed in the following 
sections of this article: Since the minimum gap behaves 
as 0(1/N), the optimal run-time in order to stay in the 
ground state scales polynomially for the Ising model. 



C. Mixed Case 

Looking at Fig. (|4]), it seems that a symmetry-breaking 
quantum phase transition typically corresponds to a 
second-order phase transition, but there are counter- 
examples: Consider more complicated energy landscape 
[lH in Fig. (j6]). In spite of the mirror symmetry of the 
energy landscape, there is a tunneling barrier throughout 
the interpolation. An analytic example for a symmetry- 
breaking first-order transition [l8| is given by a combina- 
tion of the initial Hamiltonian from the Grover problem 
with the final Hamiltonian of the Ising model 



Hn = l- |in) 



H f = 



1 N 



(7) 



where Hf has been shifted and scaled in order to preserve 
the positive definiteness. Even though this Hamiltonian 
possesses the same bit-flip symmetry as the Ising model, 
its level structure displays an avoided-level crossing at the 
critical point, i.e., it corresponds to a first-order phase 
transition with a jump between the initial and the final 
ground state(s), see Fig. (J7J. 

It can also be shown analytically that the fundamental 
gap of the combined Hamiltonian H(g) = (1 — g)Ho + 
gHf vanishes exponentially with the system size, i.e. the 
number of qubits, see, e.g., (22l [23j. 



III. DECOHERENCE IN THE ADIABATIC 
LIMIT 

In all of the above examples, we have seen that at 
the critical point, at least some energy levels become ar- 
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Figure 7: (Color Online) Sketch of the lowest eigenvalues of 
the Hamiltonian lf?Jl. One can clearly see that the spectrum 
displays an avoided-level crossing at the critical point - thus 
corresponding to a first-order transition 



bitrarily close and thus, the response times diverge (in 
the continuum limit). Consequently, during the sweep 
through such a phase transition by means of a time- 
dependent external parameter, small external perturba- 
tions or internal fluctuations become strongly amplified 
- leading to many interesting effects, see, e.g., |24l430j . 
One of them is the anomalously high susceptibility to de- 
coherence (see also [3l|): Due to the convergence of the 
energy levels at the critical point, even low-energy modes 
of the environment may cause excitations and thus per- 
turb the system. Based on the similarity between the 
quantum adiabatic algorithms and critical phenomena, 
we have argued that adiabatic quantum algorithms corre- 
sponding to the higher-order quantum phase transitions 
should be advantageous in comparison to those of first or- 
der for closed quantum systems. The present paper aims 
at generalizations to these findings when the impact of 
decoherence is considered. 

In order to study the impact of decoherence, we con- 
sider an open system described by the total Hamiltonian 
H(t) which can be split up into that of the closed system 
H sys and the bath i?bath acting on independent Hilbert 



spaces n sys <E> "H ba th = H 

H(t) = ff sys (t) + ff, 



'bath 



A H\ r 



(8) 



plus an interaction A H luti between the two, which is 
supposed to be weak A <C 1 in the sense that it does 
not perturb the state of the system drastically. Note 
that the change of the bath caused by the interac- 
tion need not be small. To describe the evolution 
of the combined quantum state |3>(t)) € H, we expand 
it using the instantaneous system energy eigenbasis 



H sys (t) |V« (*)> = E s {t) \Mt)), via 



|*(t)> = 5>»(*w»(*)> ® M*)> , 



(9) 



where a s are the corresponding amplitudes and 
\a s ) € "Hbath denote the associated (normalized but not 
necessary orthogonal) states of the reservoir. Inser- 
tion of this expansion into the Schrodinger equation 

= H(t) |$(t)) yields (h = 1) 



d_ 
dt 



(jjs\ H sys |VV 
AE sr 



r ) -t(a s \(^ s \XH int \^ r )\a r ) (10) 



with the energy gaps AE sr (t) = E s (t) — E r (t) of the sys- 
tem and the total phase (including the Berry phase) 



<p,(t) =-] dt'[E s {t') + ff&hCtO + Afl£t(*0 

-t(A(t')\Mt')) - i(us(t')\a s (t'))} (11) 



with the energy shift H?L — 



(4>s\ H int \lp r 



i?bath = ( a s \ Hbath\otr) ■ Evidently, there are two con- 
tributions for transitions in the Hilbert space % sys of the 
system: 

- The first term on the right-hand side of Eqn. (fTTJ]) 
describes the transitions caused by a non-adiabatic 
evolution [32]. Note, however, that the factor 
(a s |a r ) and the additional phases in Eqn. (ITT1) give 
rise to modifications in the adiabatic expansion. 



- The last term in Eqn. (ri0|) directly corresponds to 
transitions caused by the interaction of the quan- 
tum system with its environment. 

Since we are mainly interested in the impact of the cou- 
pling to the bath, we shall assume a perfectly adiabatic 
evolution of the system itself, i.e., without the coupling 
to the environment A = 0, the system would stay in its 
ground state. Thus, the only decoherence channel avail- 
able is heating (i.e., excitations), the phase damping and 
decay channels, for example, play no major role here. 
Considering the adiabatic condition 



(V> s |# S ys|Vv} < {AE sr y 



(12) 



the first term in Eqn. (|TU1) is negligible and the second 
one dominates. Starting in the system's ground state 
do(t = 0) = 1, which is relevant for adiabatic quantum 
computation, the excitations s > caused by the weak 
interaction AiJ; nt with the bath 



2l s = a s (T) exp{-iip s (T)} 



(13) 



can be calculated via response theory, i.e., the solution 
of Eqn. (fTU|) is to first order in A -C 1 given by 



21., w -i I dt e lA ^ (a.| (^,| Afffat |^o) K> , (14) 
h 
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where Atp rs = tp r (t) — (p s (t) . This is rather a general re- 
sult. 

In the following, after a brief review of the impact of 
decoherence on the sweep through a first-order quantum 
phase transition [33[ in section IIV1 we study in section 
IVl the impact of decoherence due to a general reservoir 
for the quantum Ising chain in a transverse field, which 
is considered a prototypical example for a second-order 
quantum phase transition. 



IV. DECOHERENCE IN THE ADIABATIC 
GROVER SEARCH 

Let us consider the Grover model ([3]) - weakly coupled 
to a bath, where we can assume the following expansion 
of the interaction Hamiltonian [33j 

N N 

XH int =\J2<7j- A, + A 2 ^ <r, • B H ■ <*j + °( x3 ) > 

i=i i,j=i 

(15) 

where A -C 1 and <Tj(t) = (crj (t) , crj (t) , a* (t)) is the vec- 
tor of Pauli matrices in the interaction picture with 
the corresponding bath operators Aj(t) and Bij(t), etc. 
Recalling the adiabatic version of Grovers search al- 
gorithm in Eqn. at the beginning of the evolu- 
tion the system has to be prepared in the ground state 
| in) = J2x=o \ x ) /V^D- This also requires that the initial 
full density operator can be initialized as a direct product 



£>(0) = gsys(O) ® £>bath(0) , 



(16) 



i.e., system and environment are not entangled at the be- 
ginning. Since in the weak-coupling limit the adiabaticity 
condition for the open system dynamics is still in leading 
order the same as for the closed system, similar to the 
discussions in the former section, one can assume perfect 
adiabatic evolution of the unperturbed system and hence 
only consider perturbations due to the interaction with 
the environment. 

The spectrum of the Grovers Hamiltonian @ con- 
sists of the ground state \ipo(g)) and the first excited 

state |V>i (<?)), which come very close \AE m i a = l/V^Dj 

at g c — 1/2, whereas all other states \ipk>i{g)) are de- 
generate and well separated from the ground state by 
an energy gap of order one. Since the temperature and 
hence the energies available in the environment must be 
much smaller than that gap of order one (in order to 
prepare the initial ground state), transitions from the 
ground state to these states \ipk>i(g)} are exponentially 
suppressed. Thus, the final probability of the transitions 
to the first excited state [33| 



N 



T dt x f dt 2 (^(ti)^(ta)) x 



(17) 



provides a good measure for the success probability which 
corresponds to |2l s | 2 <§: 1. It can be shown [33| that the 
contributions proportional to Bij do not contribute to 
second order in A. In this equation, \w) denotes the 
marked state for Grovers problem, see Eqn. ([3]), and 
is the state orthogonal to \w) in the subspace 
spanned by | w) and | in) . The expression (|17|) demon- 
strates that both system and reservoir properties affect 
the excitation amplitude. Of the system matrix elements 
(w^~\(Tj(£) | to) only those with /i = x,z contribute (for 
large D = 2 N ^> 1, the \i = y term is suppressed by a 
factor 1/y/D) 



(w x \*?{t)\w) 



i_7 g(t) 

DAE(t) 



cxp < —i 



dt'AE(t') 



(18) 



for large TV and AE(t) = y/l - 4g(t)[l - g(t)](l - 1/D). 
It is the same for (ur 1 1 cr|(t) \ w) apart from an additional 
sign (—l) Wj+1 , where Wj is the j-th bit of w, i.e., |to) is an 
eigenstate of the operators cr| with eigenvalues (—1)™^. 

Assuming a stationary reservoir [if bath) f?bath] = 
(which does not necessarily imply a bath in thermal equi- 
librium) allows for a Fourier decomposition of the bath 
correlation function 



<^(tx)^(t 2 )> = / due-***-**/ (19) 



where fj^ (cj) depends on the spectral distribution of the 
bath modes and the temperature, etc. 

For example, for a bosonic bath in thermal equilibrium 



k4 



(as e.g., 



and coupling operators oc ^ fc hkdu 

used in the spin-boson model) we would for an inverse 
bath temperature (3 obtain a Fourier decomposition such 
as 201 



J{\u\) 



1 



+ e(w) 



(20) 



where Q(lj) is the step function being 1 for uj > and 
for uj < 0. J(uS) denotes the spectral density, which is 
often parameterized as [ll| 



J(u) = Zdul^uj'e-^' 



(21) 



where < e < 1 corresponds to the sub-ohmic, e = 1 to 
the ohmic, and e > 1 to the super-ohmic case. 
Insertion of (QJ]) and ([19]) into Eqn. jT^I yields 



121., 



N 



dw fW^) 



dt 



L 9 ^ expiiut + i f dt'AEtt') 
DAE{t) I Jo 



(22) 
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plus similar terms including f^ z , f£? and ff? with the 
associated signs (—1) and (—1)^ for x and z, respectively 
[33l ] . In order to evaluate the time integrations, it is useful 
to distinguish different domains of uj: 

- For large frequencies \w\ >• Ai? m ; n , the time in- 
tegral can be calculated via the saddle-point ap- 
proximation. The saddle points t* are given by a 
vanishing derivative of the exponent 



u + AE(tZ) = 0, 



(23) 



which corresponds to energy conservation. Hence 
large positive frequencies u> 3> A_E m i n do not con- 
tribute at all which is quite natural (this corre- 
sponds to the transfer of a large energy from the 
system to the reservoir). 

The saddle-point approximation cannot be applied 
for small frequencies oj — 0(AE ul i n ) and energy 
conservation is also not well-defined. In this case, 
one might estimate an upper bound for the time 
integral by omitting all phases. 



These, altogether yield 



|a. 



X 2 D 



+AB, 
A£ mi 

2 r 1 



dujf(uj) 



duj 



/(- 



/ (MjJ — 71 — ; r , 



(24) 



where f(ui) is the appropriate sum of the ff?, f£?, 
f£? and f£? contributions. The second term of the 
above equation depends on the interpolation function 
g(t). Considering three scenarios [Hf 

(a) g = (b) gocAE (c) g cx AE 2 



the second integrand scales as 



(a) 



Df(- 



(b) 



Df(-u) 



uj- 



(c) 



/(-") 



respectively. In all of these cases, the bath modes with 
large frequencies |w| 3> AE m in do not cause problems in 
the large- N(D) limit, since the spectral function f(—oj) 
is supposed to decrease for large |w| as the bath does not 
contain excitations with large energies - the environment 
is cold enough, compare also Eqn. (|20[) . Therefore, the 
low-energy modes of the reservoir uj = 0(AE m [ n ) give the 
potentially dangerous contributions. Independent of the 
dynamics g(t) both the first integral and the lower limit 
of the second integral yield the same order of magnitude 

M 



1 3s 



A 



2 / P(AE n 



AE„ 



(25) 



as uj or even faster in order to keep the error |2t s | un- 
der control. Thus, one can conclude that the spectral 
function f(uj) of the bath provides a criterion to favor or 
disfavor certain physical implementations. If f(u)) van- 
ishes in the infrared limit faster than uj - compare also 
Eqn. (|2ip - the computational error does not grow with 
increasing system size - the quantum computer is scal- 
able. This result has already been derived in [33j with a 
slightly different formalism. 



V. RESULTS: DECOHERENCE IN THE 
TRANSVERSE ISING CHAIN 

As we shall see below, the situation may be very dif- 
ferent for second-order transitions compared to the first- 
order transitions. These investigations are particularly 
relevant in view of the announcement (see, e.g., the dis- 
cussion in [l5j]) regarding the construction of an adiabatic 
quantum computer with 16 qubits in the form of a two- 
dimensional Ising model. 

First of all, we briefly review the main steps [l| of the 
analytic diagonalization of H sys , where we switch tem- 
porarily to the Heisenberg picture for convenience: The 
set of iV qubits in Eqn. ([5]) can be mapped to a system 
of N spinless fermions Cj via the Jordan- Wigner trans- 
formation [36} given by 



0j(t) 



l-2ct(t) Cj (t), 
-II [l-2c^)Q(i) 



Cj (t) + c](t) (26) 



with (Tj(t) indicating the Pauli operators in Heisenberg 
picture (Tj{t) = Ul ys (t) <Tj £4 ys (t), where U sys {t) is the 
unitary time evolution operator of the system. It is easy 
to verify that the fermionic operators anti-commutation 
relations satisfy 



{c*,ct} 



= 5 f . 



{c t ,c j } = {4,c]}=0. (27) 



Insertion of Eqn. (|26|) into the system Hamiltonian in 
Eqn. ([5]) yields in the subspace of an even particle number 



N 



H sys (t) = - }j { [1 - 9(t)] (l - 2c] c,-) (28) 
+ 9 W (cj+i Cj + c] +1 Cj + cj Cj+i + c] c t +1 ) I , 

where the time-dependency of the Cj has been dropped 
for brevity. This fermionic Hamiltonian has terms that 
violate the fermion conservation number, Cj+\ Cj and 
Cj Cj +1 . This bilinear form can now be diagonalized by a 
Fourier transformation 



Since AE m i n decreases as 1/v-D in the large- D limit, the 
spectral function /(w) must vanish in the infrared limit 



Cj(t) 



c k (t) e 



ik(ja) 



(29) 
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followed by a Bogoliubov transformation [37j | . Here a is 
lattice spacing. The Bogoliubov transformation 



Cfe(t) = u k (t) 7fc + iv* k (t) 7l fc 



(30) 



maps the Hamiltonian into a new set of fermionic op- 
erators 7fc whose number is conserved. The same anti- 
commutation relations as in Eqn. (f27|) are also satisfied 
by 7^ and 7^ 

{jk.nl} = $kk> , iik,ik'} = {7L7!'} = 0. (31) 

Since these fermionic operators are supposed to be time- 
independent, the Bogoliubov coefficients Uk and Vk must 
satisfy the equations of motion 

i ~dT = a k(t) u k(t) + Pk(t)vk(t) , 

= -«*(*)«*(*) + , (32) 

at 

where afe = 2 — 4g(t) cos 2 (fca/2) = 2c;(t) sin(fca). For 
an adiabatic evolution | i? sys | Vv) (AE sr ) , these 
equations of motion can be solved approximately 



u k (t) « ^±^lexp^ ; / ,//'/-,.(/', }, . 



«*(*) « ex P ^ / ,//'/;,, i/'i 



(33) 



with the normalization Af k = \/2<§ k + 2otk$k ensuring 
|wfc| 2 + |«fe| 2 = 1 and the single-particle energies 

Skit) = 2y/l - 4 5 (t) [1 - g(t)} cos 2 (ka/2) . (34) 



All the excitation energies §k take their minimum 
values <?™ m = 2 |sin (ka/2)\, at the critical point 
,9cr = 1/2. The pseudo-momenta ka take half-integer 
values he (1 + 2Z)ir/N : \ka\ <ir. In view of the 
fc-spectrum the minimal gap between the ground state 
and the first excited state scales as AE m - m = 0(1/N). 
Finally, the Hamiltonian (|S|) in the subspace of an even 
number of quasi-particles reads 



lllk 



1 



(35) 



with fermionic creation and annihilation operators 7^ , jk ■ 
Hence, its (instantaneous) ground state contains no 
fermionic quasi-particles Vfc : jk IV'o(i)) = 0. Without 
the environment, the number of fermionic quasi-particles 
7feTfc wou ld be conserved and the system would stay in 
an eigenstate (e.g., ground state) for an adiabatic evo- 
lution. The coupling to the bath, however, may cause 
excitations and thus the creation of quasi-particles due 
to decoherence. 



Of course, the impact of decoherence depends on the 
properties of the bath and its interaction with the system 
(decoherence channels). In the following, we study three 
different decoherence channels. However, in all of these 
different cases, we do not specify the bath -ffbath in much 
detail for the purpose of deriving generally applicable 
results. 



A. Uniform Coupling Strengths 

Let us first consider an interaction XH m t which is 
always present: In the Hamiltonian H sys in Eqn. ([3]), 
the transverse field B(t) = 1 — g(t) appears as a classi- 
cal control parameter B c \. However, the external field 
B — > B c \ + SB does also possess (quantum) fluctuations 
SB, which couple to the system of Ising spins. Therefore, 
we start with the following interaction Hamiltonian 



SB, 



(36) 



where SB denotes the bath operator. Note that this per- 
turbation should be considered as mild, since it does not 
even destroy the bitflip symmetry of the Ising model §5§ 
and thus does not lead to leakage between the two 
subspaces of even and odd bitflip symmetry (or quasi- 
particle number, respectively). This interaction Hamil- 
tonian yields the same matrix elements as the non- 
adiabatic corrections (ip s \ H sys \ip r ) in Eqn. (|10p . which 
can therefore be calculated analogously. Insertion of 
Aifint into Eqn. (JTU) yields 

2l s « -iX [ duf s (uj) [ dt (^a(f)lX}°i M>(*)> x 



expji -ojt + J dt' AE sQ {t') 



(37) 



We may also here include all relevant properties of 
the environment into the single-operator (compare with 
Eqn. (fl9|) for the double operator version) spectral func- 
tion f(uj) of the bath 



- 4A ^<*> (a s (t)\SB(t)\a (t)) 



dioe-^Mio), 



(38) 



where Atp' s coincides with tp s — tpo in Eqn. (|11[) apart from 
the system's energy gap AE s q and is typically dominated 
by the contribution from ffbath ~ -^bath- Note that Eqn. 
(f3"5| is the generalization of Eqn. (JTHJ) for the case that 
the bath state changes strongly. As a first approxima- 
tion, we assume that f(ui) does not change significantly 
if we increase the system size N (scaling limit). After 
inserting the Jordan- Wigner, Fourier, and Bogoliubov- 
transformations, the matrix element in Eqn. (|37p reads 
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03 > 



9: 



0=1/2 



g: 



9 



a) < 



Figure 8: (Color Online) Sketch of the excitation spectrum of 
the Ising chain H aya as a function of g. For a given frequency 
oj > 0, real saddle points correspond to intersections of the 
(solid) energy level curves (e.g., AEi) with the (dashed) ver- 
tical w-line which occur shortly before (g+) and after (gt) 
the quantum phase transition at g CI — 1/2. The saddle-point 
approximation can only be applied if the intersection angle is 
large enough, i.e., for the drawn u: > line, it would work for 
AEi , but not for AE 3 etc. 



£<^K(t)|^ ) 



2ig(t) sin(fca) 



(39) 



where the ~ sign refers to the adiabatic approximation. 
Thus, it is only non- vanishing for excited states \ip s ) con- 
taining two quasi-particles s — (k, —k) with opposite mo- 
menta and hence we get AE s q = 2<§ k . 

First of all, in order to have a quantum phase tran- 
sition (or a working adiabatic quantum computer), the 
environment should be cold enough to permit the prepa- 
ration of the system in the initial ground state such that 
f(u>) is only non-negligible when 



uj < 2 = S k (t = 0) 



(40) 



holds, compare also Eqn. (|20"]) . Therefore, we will analyze 
the spectral excitation amplitude 21^ defined via 



51, 



duf s {w)W s 



(41) 



in the different w-regimes in the following. 



For the exponent in Eqn. (|37)) the saddle-point condition 
reads 



dh k (t,u) 



= ~> ijj = AE s0 {U) 
t=t, 



2S k {U), (43) 



dt 

where h k {t,uj) = i —ut + 2 JJ dt' S k (t') and denotes 

the saddle points. This condition yields two saddle points 
shortly before and after the transition, see also Fig. <j8j) 



g(tt 



l 

2* 



[uj 2 - 16 sin 2 (jfco/2)] 
8cos(fca/2) 



1/2 



(44) 



The saddle-point approximation yields for the spectral 
excitation amplitude defined in (I4ip 



2^oj»2|fea| 



±32iriX 2 sin(fea) sin (if) e 2 M*.,") 
\\ ug{U)g-\U)^rf -16 sin 2 (if) 



-O 



\g(t* 



U! V U! 2 



Ak 2 o? 



(45) 



which depends on the interpolation dynamics g(t). The 
minimum gap can be obtained from Eqn. (1341) and does 
indeed scale polynomially AE nlin — 0(1/N) and, thus: 

- For a constant speed interpolation g(t) = t/T, the 
necessary run-time for an adiabatic evolution T 
scales polynomially T = O (AE~ 2 n ) = 0{N 2 ). 

- For adapted interpolation dynamics g(t) cx AE(t) 
or g(t) cx AE 2 (t), however, one may achieve shorter 
run-times of T = O(NlnN) or T = O(N), respec- 
tively [34| and therefore better results for the spec- 
tral excitation amplitude, see Table |U 



2. Near the Minimum Gap 

From Eqn. ([4"5"j) it follows that the saddle-point approx- 
imation breaks down if w approaches the minimum gap 
w s» AE^ m 2\ka\, see Fig. ©. In this case, we may 
obtain an upper bound for the time integral in Eqn. (|37[) 
via omitting all phases. For a constant speed interpola- 
tion g(t) = t/T 



Sit 



= 2|/vrj 



< 



2X sin(fca) 
T 

O (\N 2 uj\nuj) 



t 



T 

dt- 

(i <Sk\t) 



(46) 



Similarly, one can get better results for adapted interpo- 
lation dynamics, see Table HI 



1. Intermediate Positive Frequencies 

We may solve the time integral via the saddle-point (or 
stationary phase) approximation for intermediate posi- 
tive frequencies, 



2 » cj » AE™ n 



2\ka\ 



(42) 



3. Positive Frequencies Below the Minimum Gap 

For positive frequencies fulfilling < u <C 2|fca|, the 
saddle points at 



g(tt 



l . l 



± -\/uj 2 - Ak 2 a 2 , 



(47) 
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1 > W > 2fca 


1 > OJ « 2fca 


g(t) = o 

g(t) oc AE(t) 
S^cx A£ 2 (i) 


O (Xkauj- 2 ) 


C(AiV 2 tjlntj) 
O(AiVlnjV) 
O(AiV) 



Table I: Scaling of the spectral excitation amplitude 21" in the 
saddle-point approximation (uj ^> 2ka) and its upper bound 
(id ~ 2fca) for different interpolation dynamics g(t), where 
AE(t) — 2<f fc= ^/( a jv) (t) denotes the fundamental gap. In all 
cases, the total excitation probability (integral over all ui and 
sum over all k) increases with system size N. 



move away from the real axis and thus the exponent in 
Eqn. p7p contains real terms. The constant speed inter- 
polation leads to 



2 / dtS k (t) 
'o 



IT] 



uj\ka\ (kaY 



T. 



Im(t) 




(48) 



Figure 9: (Color Online) A sketch of the deformed integration 
contour. The original integration contour (blue line along the 
real axis) is shifted to the complex plane (curved line), where 
t, indicates the singular point, t, = T/2 + iT/2 tan (ka/2). 
Only paths a and b contribute significantly to the integral. 



where rj is a real value. Therefore, the spectral excitation 
amplitude is exponentially suppressed in the adiabatic 
limit 



Of 



cj<§C2| ka\ 



= 0[r y . v {~T(ka) 



(49) 



4- Negative Frequencies 



Finally, for negative frequencies uj < 0, the saddle 
points collide with the branch cut generated by the 
square-root in S k . In this case, we may also estimate the 
spectral excitation amplitude 21" in Eqn. (14TT) by deform- 
ing the time integration contour into the complex plane. 
We assume that all involved functions can be analytically 
continued into the complex plane and are well-behaved 
near the real axis. Given this assumption, we deform the 
integration contour into the upper complex half-plane to 
obtain a negative exponent which is the usual procedure 
in such estimates until reaching a saddle point, a singu- 
larity, or a brunch cut, see Fig. ©. Deforming the inte- 
gration contour into the lower complex half-plane would 
of course not change the result, but there the integrand 
is exponentially large and strongly oscillating such that 
the integral is hard to estimate. Since the integral in the 
complex plane is zero around path c and the integrals on 
the paths 1 and 2 cancel each other, only paths a and b 
give the main contribution to the integral. 

Let us first consider a constant interpolation function 
g(t) = t/T which leads to singular points 



T T fka\ 
t* = — ±i— tan — 
2 2 V 2 J 



(50) 



in the complex plane. Performing the time integral in 
the exponent of Eqn. (|3"T|) acquires a large negative real 



term in the exponent 



dt' AE s0 (t') 




£'±2t 2 + ^(fca) 2 , 



dt' AE s0 (t') 



(51) 



where £' is a constant and real value. Insertion of 
Eqn. ([STjl into Eqn. (|3T|) and doing some algebra yields 
the exponential suppression for the amplitudes in the up- 
per complex half-plane 



51 



u:<0 



dt 



a . b 



F(t) exp ji ( 



2r -ut + T 



) } (52) 



with F(t) — (ip s (t)\J2j Tj \ipo(t)} and where r is a real 
constant. Therefore, applying the inequality 

\Sj\< f dy\<f>(y)\ with Sj = J dy <j>{y) , (53) 

the amplitudes for negative frequencies are also expo- 
nentially suppressed for g(t) = t/T and similarly for the 
other interpolations. This result can be understood in the 
following way: For frequencies uj below the lowest excita- 
tion energies, the energy lj of the reservoir modes is not 
sufficient for exciting the system via energy-conserving 
transitions. Hence excitations can only occur via non- 
adiabatic processes for which energy-conservation be- 
comes ill-defined, but these processes are suppressed if 
the evolution is slow enough. 

An estimate on the total error probability |2l s | 2 in- 
troduced in Eqn. (I41[) is obtained by performing the 
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weighted sum of the contributions from the different ui- 
regimes, which depends on the Fourier transform of the 
bath correlation function 



I at 



< 



max 

uj^>2\ka\ 



21 



lo^>2\ ka\ 



\f s {uj)\ dui 



2»u»2|/ca| 



max 

jK2\ka\ 



^uj~2\ka\ 



|/ s (w)| du) 



2\ka\ 



max 

0<ui<2\ka\ 



2pj<2|fca| 



+ 



max|2l^ <0 | 



\f s {u)\dw . 



(54) 



From the results in Eqns. (g5]), gS]), (gHJ), and ^ it be- 
comes obvious that although the last two contributions 
in the above sum can be efficiently suppressed, the first 
two terms will scale with the system size. Therefore, for 
the adiabatic Ising model, decoherence can only be effec- 
tively suppressed when the bath spectral function / s (cj) 
has only support at frequencies below the minimum fun- 
damental energy gap. For a bosonic bath in thermal 
equilibrium this would imply a reservoir temperature be- 
low the minimum fundamental energy gap, compare also 
Eqn. J20]). 



B. Nonuniform Coupling Strengths 

In a realistic situation, the Ising spins will not be sym- 
metrically coupled to the environment 



H in t = E ® SB j , 



(55) 



where 5Bj r = (a s (t) \ SBj(t) \a r (t)), A(p' s coincides with 
ip s — ipo in Eqn. (TTTT) apart from the system's energy 
gap AE s0 and u k (t) , v k (t) are the Bogoliubov coefficients 
given in Eqn. (|33[) . The excitation amplitude in Eqn. (|56[) 
is only non-vanishing for the excited states containing 
two quasi-particles s = (fc, — k') 



(58) 



A£ s0 = S k + ■ 



Insertion of u k , v k ' given by Eqn. (|33p . into Eqn. ([5 
yields 



a. 



k,k' 



exp 



-uit - 



2 / dt' ' Sk(t') 
/o 



(59) 



where 

C k ,k'(t) =4 5 (t) sin(fc'a) 




2g(t) cos 2 (fca/2) 



<&(*) 



(60) 



Following the same procedure outlined above, we may 
consider different domains of u> in order to solve the time 
integral in Eqn. (|59[) . For the intermediate positive fre- 
quencies, 

2 > ui > 2 \ka\ and 2 \k'a\ , (61) 

the saddle-point approximation can be applied once again 

-w+ 2**(i„) = 0, (62) 

where t* denotes the saddle points. The saddle-point ap- 
proximation yields for the spectral excitation amplitude 



where SBj denote now different operators acting on the 
bath Hilbert space. Insertion of (|55|) into Eqn. (fT4"f and 
evaluating the corresponding matrix element by apply- 
ing the Jordan- Wigner, Fourier, and Bogoliubov trans- 
formations of Pauli matrices - given in Eqns. (|2"6")l . (|2"9"1) . 
and (jnnj) - yield 



k.k' J Ja 



-uit 



dt u* k (t)vk'(t) x 



dt' AE s0 (t') 



(56) 



where we include again all relevant properties of the en- 
vironment into the spectral function f k k ,(ui) of the bath 



i{k—k')ja 



+oo 



dwe-^f^iu), 

(57) 



93 



lo^>2\ ka\ 



o 



Xk'a 



\J ui g(t*)\fu) 2 — Ak 2 a 2 
where 03^ is defined as following 



(63) 



21., 



E 

k,k' 



du ft 



(64) 



The spectral excitation amplitude in Eqn. (|63|) depends 
on the interpolation dynamics g(t) and is very similar to 
Eqn. (|45j) apart from N in the denominator. Presence 
of N in the denominator may cause some slow down for 
the error probability, see Table HI However, existence 
of many excited states - sum over all possible k and 
k' in Eqn. (|59"|) - causes the growth of the error prob- 
ability with the system size. If ui approaches the 2|fca|, 
the saddle-point approximation breaks down and we can 
get a similar upper bound which is shown in Table. U 
by omitting all phases. With the same argument given 
above, the amplitudes are exponentially suppressed in 
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the adiabatic limit for frequencies far below 2|fca|. Thus, 
as with the discussion in the previous subsection below 
Eqn. (|54j) we may conclude that in a more general case 
where the coupling strength to the bath is not uniform, 
the impact of decoherence for the environmental noise is 
very similar to the coherent reservoir and the error prob- 
ability increases with the system size, unless the bath 
temperature lies below the minimum energy gap. 



C. Perturbing the Bitflip Symmetry 

Unfortunately, interactions with the reservoir cannot 
be tailored, such that one may also expect perturba- 
tions that do not reflect the bitflip symmetry of the Ising 
Hamiltonian and thereby lead to leakage between the 
subspaces of even and odd quasi-particle numbers. Let 
us consider a simple case, where only one Ising spin is 
coupled to the environment 



(65) 



Hi nt = al®5B. 



Insertion of this interaction Hamiltonian into Eqn. (| 14[) 
and evaluating the corresponding matrix element yield 
the excitation amplitude which is a combination of two 
terms 



a. 



iX 



y^/.MX^ + ajk) , (66) 



with 
and 



•5,1 



•A'fc 



(67) 







g(t) = o 

g(t) oc AE(t) 
g(t) « AE 2 (t) 


O (xka/V^tj H 


+ Oi 

-o(x 
+ o 


AiV 3 / 2 ) 
vWlniv) 
(Av^v) 



Table II: Scaling of the spectral excitation amplitude 
21" for different interpolation dynamics g(t), where 
AE(t) = (ok=Tr/(aN)(t) denotes the fundamental gap. In all 
cases, the total excitation probability (integral over all ui and 
sum over all k) increases with system size N. 



= f_ T dta(t) e~ luJt , and f2(w) is some function of 
lj. In order to solve the time integral of 21" 2 , we may 
employ saddle point approximation which yields 



Thus, we can conclude that 



2L 



du> f s (u) 



E 

k ■ 



(71) 



(72) 



The optimal run-time for an adiabatic evolution T de- 
pends on the interpolation dynamics g(t), see 34|. Scal- 
ing of the spectral excitation amplitude 21" for differ- 
ent interpolation dynamics is shown in Table HU This 
decoherence channel poses a significant problem to ro- 
bust ground state preparation in the Ising model: Re- 
gardless how low the bath temperature is, the decoher- 
ent excitation probability will scale with the system size! 
This finding is consistent with the experience that large 
Schrodinger cat states as ^ are extremely sensitive to 
decoherence. 



2lt 



1 l-2 g ft)cosW2) x 
V 2 &(t) 



x (V ) s (*)l7fclV'o(t)>expii 



-ujt + 2 / dt' ' £ k {t') 
Jo 



(68) 



where the spectral function f a (w) of the bath is as defined 
in Eqn. (|55| . For large T, the first term 21" 1 is a Fourier 
transformation of some function of uj 



( 



E(w) - I dt H(t)e" 

■T 



O(l/o>) 



t 

J 



where 

m = 



y/2££(t) + 4 [1 - 2g(t) cos 2 (fca/2)] S k (t) 



(69) 
, (70) 



VI. SUMMARY 

In summary, we studied the impact of decoherence due 
to a weak coupling to a rather general environment on 
first and second order quantum phase transitions. Since 
the Ising model is considered [1| as a prototypical ex- 
ample for a second-order quantum phase transition, we 
expect our results to reflect general features of second- 
order transitions. 

For the decoherence channel (|36[) which is always 
present (though possibly not the dominant channel), we 
already found that the total excitation probability al- 
ways increases with system size ./V (continuum limit): 
Even though the probability for the lowest excitation 
k = ±ir/(aN) can be kept under control for a bath which 
is well-behaved in the infra-red limit - see also Sec. (3.4) - 
the existence of many excited states k £ n(l+2Z) / (aN) : 
\ka\ < 7r converging near the critical point causes the 
growth of the error probability for large systems. This 
growth can be slowed down a bit via adapted interpola- 
tion schemes <?(t), but not stopped. 
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Other decoherence channels will in the best case dis- 
play the same general behavior: E.g., for 2 3> ui 3> \ka\, 
the associated amplitudes scale as 



Xcj) s (u 



(73) 



where <p s denotes the matrix element in analogy to 
Typically, for a homogeneous coupling to the bath, <j> s 
does not strongly depend on the system size N (for given 
ka and U)). Since <?(i*) decreases for N — > oo or at least 
remains constant - for g(t) oc AE 2 (t) - the total excita- 
tion probability again increases with system size N. If 
only a few spins are coupled via crj to the environment, 
the matrix element (f> s in Eqn. (j5"o| will decrease with the 
system size 0(1/N) and then the error probability may 
be kept under control - for g(t) oc AE 2 (t). However, 
4> s is of order 0(1/ y/N) for the the cr z -channel given in 
Eqn. (|66l) and the Schrodinger cat states are still sensitive 
to decoherence. 

According to the analogy between adiabatic quantum 
algorithms and quantum phase transitions @, Q , this re- 
sult suggests scalability problems of the corresponding 
adiabatic quantum algorithm - unless the temperature 
of the bath stays below the ( A-dependent) minimum gap 
[H or the coupling to the bath decreases with increasing 
A. These problems are caused by the accumulation of 
many levels at the critical point g cr = 1/2, which presents 
the main difference to isolated avoided level crossings 
(corresponding to first-order phase transitions) discussed 
earlier. It also causes some difficulties for the idea of 
thermally assisted quantum computation (see, e.g., [H 



since, in the presence of too many available levels, the 
probability of hitting the ground state becomes small. 

Therefore, in order to construct a scalable adiabatic 
quantum algorithm in analogy to the Ising model, suit- 
able error-correction methods will be required. As one 
possibility, one might exploit the quantum Zeno effect 
and suppress transitions in the system by constantly mea- 
suring the energy, see for example [39| . Another interest- 
ing idea are adiabatic ground state preparation schemes 
(algorithms) that provide a constant lower bound on the 
fundamental energy gap that does not scale with the 
system size. In case of the Ising model discussed here, 
this is possible by increasing the complexity of the in- 
terpolation path (i.e., beyond the straight- line interpola- 
tion). Unfortunately, the simplest approach to the Ising 
model [401 ] only bounds the fundamental gap in the sub- 
space of even bitflip parity, i.e., decoherence channels 
that mediate transitions between the two subspaces as 
e.g. in Eqn. (|65|) will destroy the associated robustness 
against decoherence. Many-body interactions in the sys- 
tem Hamiltonian are required to resolve this problem. In 
this case, decoherence could be strongly suppressed for 
a low-temperature bath. Of course, the generalization of 
all these concepts and results to more interesting cases 
such as the (NP-complete) two-dimensional Ising model 
is highly non-trivial and requires further investigations. 
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